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Compton Deconvolution Camera 
[0001] This application claims the benefit under 35 U.S.C. § 1 19(e) of U.S. Provisional 
Patent application no. 60/205,796, filed on May 17, 2000, and entitled, "Compton 
Deconvolution Camera". 
FIELD OF THE INVENTION 

[0002] The present invention pertains to medical imaging systems. More particularly, 
the present invention relates to a method and apparatus for generating images in a 
medical imaging system based on detected radiation by applying Compton scattering 
and deconvolution principles. 
BACKGROUND OF THE INVENTION 

[0003] In nuclear medicine, images of internal structures or functions of a patient's 
body are acquired by using one or more gamma camera detector heads to detect 
radiation emitted by a radiopharmaceutical that has been injected into the body. A 
computer system controls the gamma camera detector heads to acquire data and then 
processes the acquired data to generate images. Nuclear medicine imaging techniques 
include single photon emission computed tomography (SPECT) and positron emission 
tomography (PET). SPECT imaging is based on the detection of individual gamma rays 
emitted from the body, while PET imaging is based on the detection of gamma ray pairs 
resulting from electron-positron annihilations and emitted in coincidence with each 
other. 
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[0004] Existing gamma cameras commonly utilize a heavy lead collimator, a planar 
scintillation crystal, and a coarse array of bulky photomultiplier tubes (PMTs). The 
collimator absorbs all incident gamma photons except those traveling in a specific 
direction at each cell in the collimator in order to reconstruct a meaningful two- 
dimensional planar image. Collimators with different acceptance geometries are 
available, but must be physically exchanged with a cumbersome transport mechanism 
capable of supporting their substantial weight. 

[00051 A collimator drastically reduces the counting efficiency of the gamma camera. 
It adds great weight to the camera detector heads which, in turn, must be supported by 
a very heavy gantry. Exchanging collimators requires time and labor, thereby reducing 
patient throughput. The bulky exchange mechanism occupies a sizeable space and is 
subject to mechanical misalignment, thereby requiring service calls. 

[0006] PMTs are an old technology that has not fundamentally changed in decades. 
PMTs are not only bulky, but also consume energy, generate heat, and require 
replacement when their efficiency decays or when they fail. 
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SUMMARY OF THE INVENTION 

[0007] The present invention includes a Compton deconvolution camera that 
comprises at least two detection layers, each to detect events resulting from incident 
photons. The camera further includes position sensing logic to determine positions of 
events in each of the detection layers, a coincidence detector to detect pairs of coincident 
events resulting from Compton scattering, and processing logic. For each of a plurality 
of subsets of one of the detection layers, the processing logic associates data 
representing detected events with a distribution of corresponding events in the other 
detection layer, based on the detected pairs of coincident events. The processing logic 
then uses a deconvolution function to localize probable source locations of incident 
photons, and then uses the probable source locations to reconstruct an image of an 
object. 

[0008] Other features of the present invention will be apparent from the 
accompanying drawings and from the detailed description which follows. 
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BRIEF DESCRIPTION OF THE DRAWINGS 

[0009] The present invention is illustrated by way of example and not limitation in the 
figures of the accompanying drawings, in which like references indicate similar 
elements and in which: 

[0010] Figure 1 A illustrates a biplane scintillation camera with arrays of solid-state 
wafer photodetectors; 

[0011] Figure IB illustrates a hybrid biplane camera with a first solid-state ionization 

detection layer and a second scintillation layer with PMTs; 

[0012] Figure 1C illustrates a biplane solid-state ionization detector; 

[0013] Figure 2 is a block diagram of a biplane coincidence gamma camera; 
[0014] Figure 3A is a side view of a biplane scintillation camera with arrays of solid- 
state wafer photodetectors; 

[0015] Figure 3B is a side view of a hybrid biplane camera with a first solid-state 
ionization detection layer and a second scintillation layer with PMTs; 

[0016] Figure 3C is a side view of a biplane solid-state ionization detector; 

[0017] Figure 4A illustrates coincidence detection of a single incident photon through 
a Compton scattering event in the first detection layer and an absorption event in the 
second detection layer, showing also the probability cone of origin; 
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[0018] Figure 4B illustrates multiple coincident events arising from the same point 
source of incident photons, whose cones-of-origin can be deconvolved to the original 
point source; 

[0019] Figure 4C illustrates the probability distribution of absorption events in the 
second detection layer resulting from incident photons from a point source undergoing 
Compton scattering at one location in the first detection layer; 

[0020] Figure 4D illustrates the combined probability distribution of absorption 
events in the second detection layer resulting from incident photons from two point 
sources undergoing Compton scattering at the same location in the first detection layer; 

[0021] Figure 4E illustrates the skewed probability distribution of absorption events 
in the second detection layer for photons incident at a non-perpendicular angle on the 
first detection layer, and the reconstruction of an unskewed distribution in a virtual 
plane perpendicular to the trajectory of the incident photons; 

[0022] Figure 5 illustrates a method of detecting multiple coincident events 
simultaneously by limiting the definition of coincidence to be within cones-of- 
acceptance; 

[0023] Figure 6A illustrates Compton scattering of incident photons within a range of 
angles about the original angle of incidence; 
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[0024] Figure 6B illustrates an angular probability distribution of detected Compton- 
scattered photons, in which failure to detect incident photons that do not scatter results 
in a hole around 0 degrees; 

[0025] Figures 6C, 6D and 6E illustrate, respectively, the inverse deconvolution 
kernel, deconvolution kernel, and resulting delta function defining the probable angular 
origin of incident photons; 

[0026] Figure 6F illustrates the calculation of angle and offset values for photon 
incidence at oblique angles on detector planes rotating about an isocenter; 

[0027] Figure 7 is a schematic drawing of the arrangement of image memory for 
associating coincident counts in the second detection layer to counts at each pixel in the 
first detection layer; 

[0028] Figure 8 is a flowchart of the operations in performing event detection and 
image reconstruction using the Compton Deconvolution Camera (CDC); and 

[0029] Figures 9A through 9F illustrate several image vector formats possible with 
variable electronic collimation. 
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DETAILED DESCRIPTION 

[0030] A method and apparatus for generating images in a medical imaging system 
based on detected radiation by using Compton scattering principles are described. An 
embodiment of the apparatus is referred to herein as the Compton Deconvolution 
Camera (CDC). Note that in this description, references to "one embodiment" or "an 
embodiment" mean that the feature being referred to is included in at least one 
embodiment of the present invention. Further, separate references to "one 
embodiment" in this description do not necessarily refer to the same embodiment; 
however, neither are such embodiments mutually exclusive, unless so stated and except 
as will be readily apparent to those skilled in the art. Thus, the present invention can 
include any variety of combinations and/or integrations of the embodiments described 
herein. 

[0031] The CDC allows imaging of the spatial distribution of radioactive sources (e.g., 
the biodistribution of radiopharmaceuticals) without a physical collimator. By 
"physical" collimator, what is meant is a piece of equipment that absorbs photons not 
traveling in a specific direction to prevent them from reaching a detector, an example of 
which is a conventional lead collimator. Although described in greater detail below, the 
CDC may be briefly summarized as follows. A biplane-detector solid-state gamma 
camera head determines the angular origin of each incoming gamma photon by 
coincident events (scintillation or ionization) in two parallel detector planes. Spatial 
localization is achieved by deconvolution of a large number of such coincident events. 
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Thus, incident gamma photons are accepted from a wide and variable field-of-view 
instead of through a collimator. 

[0032] Each detection layer includes either solid-state ionization detectors or a 
scintillation crystal plane with solid-state photodetectors. Alternatively, the second 
layer (i.e., the layer farthest from the object to be imaged) may instead have 
photomultiplier tubes (PMTs). The first layer (the layer closest to the object to be 
imaged) is optimized to increase the probability of occurrence of a single Compton 
scattering. The second layer is designed to maximize the probability of absorption of all 
residual energy in the scattered gamma photon for accurate determination of photon 
energy. 

[0033] Most image reconstruction can be done in software from coincident counts 
collected in the two detection layers. As a result, many imaging parameters can be 
adjusted in post-processing—including collimation (if desired: parallel, convergent, 
divergent), steering, angle of acceptance, coincidence filtering, and /or focusing— 
providing enormous flexibility in image reconstruction and image reprocessing. 

[0034] By using two detection planes, the angular origin of each incident gamma 
photon can be determined, eliminating the need for a collimator. The use of solid-state 
devices is independent of the collimator, but provides a way of eliminating 
inhomogeneous electron densities (such as metal electrodes in PMTs) that may 
otherwise scatter photons in between the two detection layers. 
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[0035] The use of solid-state devices offers many new opportunities, including higher 
resolution and improved photon spectrometry. Solid-state detector arrays available 
today offer practical alternatives to PMTs while providing improved energy resolution. 

[0036] The combination of no physical collimator, two detection layers, and solid- 
state devices offers simultaneous additional opportunities for unprecedented flexibility 
in image reconstruction, increasing sensitivity, reducing size and weight, improving 
patient throughput, and reducing mechanical cost. 

[0037] In one embodiment, the CDC includes two parallel detection layers, circuitry 
and algorithms to localize ionization or scintillation events in each detection layer and 
estimate the amount of energy deposited by an incident gamma photon, algorithms to 
compute the cone of origin of a Compton-scattered gamma photon, a deconvolution 
processor to localize the probable origin of each such gamma photon, and a device to 
visually display the reconstructed image. A nuclear medicine imaging system may 
contain one or more CDC detector heads, each including two detection layers and 
associated circuitry as described above, which may be rotated about a patient by a 
gantry to acquire data for tomographic imaging. 

[0038] The purpose of the two parallel detection layers is to not only detect an 
incident gamma photon and localize its position of incidence, but to determine from 
which direction the photon is coming. The first layer is designed to produce a high 
probability of detecting a single Compton scattering event, and the second layer is 
designed to produce a high probability of absorbing the residual energy of gamma 
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photons scattered from the first layer. In general, this is done by making the first layer 
thinner and less electron-dense, and by making the second layer thicker and more 
electron-dense. 

[0039] Three possible embodiments of the double-layer detector include: 1) two 
parallel scintillation layers, each with an array of solid-state photodetectors, as shown in 
Figures 1 A and 3A; 2) two parallel two-dimensional arrays of solid-state ionization 
detectors, as shown in Figures 1C and 3C; and 3) a hybrid camera which has a first 
semiconductor ionization detector and a second scintillation layer with an array of 
photodetectors, e.g. PMTs, as shown in Figures IB and 3B. Thus, Figure 1A illustrates a 
biplane scintillation camera head 100 comprising two detector layers including arrays 
126 and 136 of solid-state wafer photodetectors. As shown in Figure 1A, coupled to the 
photodetector arrays 126 and 136 are, respectively, diffraction layers 124 and 134, 
scintillation crystals 122 and 132, and reflective layers 120 and 130. A gamma photon 
110 is also shown, which causes a Compton scintillation event 112 in the first 
scintillation layer 122, which is scattered along path 114 to result in a total absorption 
scintillation event 116 in the second scintillation layer 132. Figure 3A shows a side view 
of the detector head 100. Visible in this view are the anterior electrical connection layers 
320 and 330 and the posterior electrical connection layers 324 and 334 for the solid-state 
photodetector elements 322 and 332, respectively. 

[0040] Figure IB illustrates a hybrid biplane camera head 101 with a first solid-state 
ionization detection layer and a second scintillation layer with FMTs. In Figure IB, the 
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first detector layer includes a solid-state ionization detector array 150 including a 
number of semiconductor ionization detector elements 154 and a posterior electrical 
connection layer 152. The second layer includes a scintillation crystal 132 optically 
coupled to an array of PMTs 138 via diffraction layer 134, and optically coupled to a 
reflective layer 130. A gamma photon 110 is also shown, which causes a Compton 
ionization event 112 in the first detector layer 150, followed by a total absorption 
scintillation event 116 in the second detector layer 132. Figure 3B shows a side view of 
the detector head 101. Visible in this view are the posterior electrical contacts 352 for 
the individual ionization detector elements 154, separated by spaces 350, as well as the 
anterior electrical connection layer 356 for the first detector layer. 

[0041] Figure 1C illustrates a biplane solid-state ionization detector. In this 
embodiment, the first detection layer is the same as in Figures IB and 3B. However, the 
second detection layer includes a solid-state ionization detector array 140 and a 
posterior electrical connection layer 142. Figure 3C shows a side view of the detector 
head 102. Visible in this view are the posterior electrical contacts 342 for the individual 
ionization detector elements 344 of the second detection layer, separated by spaces 350, 
as well as the anterior electrical connection layer 346 for the first detector layer. 

[0042] Figure 2 is a block diagram of a CDC system according to one embodiment. As 
shown, the CDC system includes two detection layers 212 and 214, The CDC system 
further includes (associated with the first detection layer 212) preamplification and 
analog-to-digital (A/D) conversion circuitry 222, position sensing logic 232, pulse 
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height analyzer 242, and image memory 252. In addition, the CDC system includes 
(associated with the second detection layer 214) pre-amplification and A/D conversion 
circuitry 224, position sensing logic 234, pulse height analyzer 244, and image memory 
254. A coincidence detector 236 is coupled to position sensing logic 232 and position 
sensing logic 254. Gating unit 256 is coupled to image memories 264 and 266 to 
facilitate dynamic cardiac and breathing scans. The acquisition computer 262 processes 
the data stored in memories 252 and 254 to reconstruct images of the radiation source 
202. The acquisition computer 262 may also control various other functions of the CDC 
system, including movement of the gantry (not shown) to rotate the detector heads 
around the patient. Deconvolution logic 264 is used to deconvolve data stored in image 
memory 254 (for the second detector) for each pixel of memory 252, to localize probable 
source locations of incident photons 204, Reconstructed images are stored in image 
reconstruction memory 266 and displayed on display device 272. 

[0043] In a typical embodiment, at least the detection layers 212 and 214 are 
implemented within a CDC camera " detector head" (not shown). The acquisition 
computer 262, deconvolution logic 264, and image reconstruction memory 266 are 
implemented within a conventional computer system, which may be a personal 
computer (PC) or a workstation, for example. The other components in Figure 2 may be 
distributed in essentially any manner between the detector head and the 
aforementioned computer system, or they may be implemented in one or more 
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physically separate modules (while maintaining the described functional 
"connections"). 

[0044] Note that many of the elements shown in Figure 2 may be embodied in either 
hardware or in software, or in a combination of hardware and software, unless 
otherwise stated herein, and except as will be recognized by those skilled in the relevant 
art. 

[0045] Hence, the weak signals of the solid-state detectors 212 and 214 are 
individually preamplified and A/D converted in the camera head. The channels from 
each detection layer are connected to position logic 232 or 234 (as appropriate), which 
performs two-dimensional interpolation of each ionization or scintillation event to 
estimate the center of its scintillation or ionization distribution within each detection 
layer, as shown in Figures 2 and 8. A coincidence detector (e.g., an AND gate) 
determines whether each ionization event has a twin that occurred in the other 
detection layer within a specified time (typically from 0.5 to 1 nanoseconds, depending 
on the separation between the two detection layers). Non-coincident events are 
discarded. Coincident events may also be screened on the basis of lateral proximity to 
accept only a certain range of probable acceptance/ scattering angles or to eliminate 
unrelated coincident events, as shown in Figure 5. Figure 5 illustrates detection of 
multiple coincident events simultaneously at locations 512 and 522 by limiting the 
definition of coincidence to be within cones of acceptance of coincident events. The 
incident paths 510 and 520 are shown, as are the paths 514 and 524 of the Compton- 
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scattered photons, which are absorbed in the second detection layer at points 516 and 
526 respectively. These paths lie within the respective non-overlapping cones of 
acceptance 518 and 528 of coincident events, and hence, unambiguously define two 
simultaneous coincident events. 

[0046] The energy absorbed in each ionization event is estimated by the pulse height 
analyzers 242 and 244 from spatially-integrated pulse height and known calibration 
coefficients for scintillation efficiency, ionization signal amplitude, etc. As a validity 
check on these values, the estimates of absorbed energy in the two detection layers are 
summed to determine if they match the expected energy of the incident gamma photon 
from the radionuclide being imaged, as shown in Figure 4A. In one embodiment, the 
validity check is performed by the acquisition computer 262. Figure 4A illustrates 
coincidence detection of a single incident photon through a Compton scattering event at 
location 112 in the first detection layer 122 and an absorption event at location 116 in the 
second detection layer 132, showing also the probability cone of origin 420. The cone of 
origin 420 is defined by the angle 6 about the central axis 418 from the location of event 
112. If the image acquisition is gated, only coincident events occurring within periodic 
intervals are taken. The selected events are stored as counts in memories for the 
respective detection layers. Note that the "memories" for the two detection layers may 
actually be embodied in a single memory, i.e., in the same physical storage device. 

[0047] For each coincident event pair of events, the scattering angle is computed, as 
described below. Considering all counts at a given pixel location in the first detection 
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layer for which there are acceptable coincident counts in the second detection layer, the 
positional distribution of all such coincident counts in the second layer is accumulated 
in image memory for each pixel location in the first layer, as shown in Figure 7. Figure 
7 is a schematic drawing of the arrangement of image memory for associating 
coincident counts in the second detection layer to counts at each pixel in the first 
detection layer. Specifically, the image memory 720 for the second detection layer is 
shown, which includes various planes 722 for forward-projected scattering distributions 
coincident with events recorded in a particular pixel 712 of image memory 710 of the 
first detection layer. The numerals inside the various pixels 724 of the second layer 
image memory 720 indicate an example of the relative number of counts in the second 
layer corresponding to one pixel 712 in the first layer image memory 710. 

[0048] After accumulating a large collection of coincident events, a deconvolution is 
performed on the data associated with the second detection layer for each pixel in the 
first detection layer to localize the probable origin distribution of the collection of 
incident gamma photons. Merging the probable origin data for all pixels in the first 
detection layer further improves the count statistics to improve the signal-to-noise ratio 
and resolution. 

[0049] At this point, additional filtering may be performed, which may be 
implemented in software for example, to achieve collimation, angle of acceptance, 
steering, focusing, etc., as shown in Figures 9A through 9F. Figures 9A through 9F 
illustrate several image vector formats possible with variable electronic collimation, i.e., 
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respectively: parallel rectangular (Figure 9 A), multi-directional within sector (Figure 
9B), divergent sector (Figure 9C), convergent sector (focused; Figure 9D), parallel 
steered (Figure 9E), and compound parallel steered (Figure 9F). 

[0050] Rotating the detector heads (or acquiring data from a ring of CDC detectors) 
begins the tomographic reconstruction process to further improve estimates of depth in 
three-dimensional image reconstruction. The result is plotted as counts in either two- 
dimensional planar projections or three-dimensional volumetric datasets, and displayed 
to the operator. 

[0051] The principle of operation of the CDC will now be further described. 
Radionuclides used in medical radiopharmaceuticals typically emit mono-energetic 
gamma photons in the energy range of 80 keV to 550 keV. In this range the dominant 
interactions in matter are photoelectric absorption and Compton scattering. If 
photoelectric absorption occurs in the first detection layer, there is no coincidence and 
the gamma photon is simply not counted. 

[0052] However, if a Compton scattering occurs, the incident gamma photon y 0 is 
scattered by collision with an electron, transferring a portion of its energy to the recoil 
electron. The trajectory of the scattered photon is at a specific angle 6 relative to that 
of the incident photon y 0 , as given by equation (1), where s 0 is the energy of the incident 
photon, 8^ is the energy of the scattered photon, m re is the rest mass of the recoil 
electron, and c is the velocity of light. 
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0 = cos _I [(l "-)-*— +1] (!) 

[0053] Thus, given the energies of the incident and scattered photons, the Compton 
scattering angle can be determined. The energy e re of the recoil electron, estimated by 
pulse height analysis in the first detection layer, is also constrained by equation (2): 

e re = (2) 
[0054] Therefore, the energy e re can be used as used as a validation check that two 

coincident events are in fact related to the same incident gamma photon. If this 

validation fails, the coincident events are discarded. This feature is useful in discarding 

those gamma photons undergoing multiple Compton scattering events within the first 

detection layer, whose 9 is not well defined. 

[0055] A coincident event in the two detection layers is accepted when the two events 
are separated in time by the interval x 2 - x v as given by equation (3), where x x and x 2 are 
the times of ionization events in the first and second detection layers, respectively, d 12 is 
the distance between the two detection layers, and v K is the angle of incidence of the 
scattered photon y x on the second detection layer: 

x 2 ~x 1 = d u I (c sin xj/^) (3) 
[0056] For example, for a detection layer separation of 15 cm and a minimum angle of 

incidence of 60° for the scattered photon y x , the coincidence interval will be < 0.58 nsec 

(see Figure 5). Even if the coincidence interval is taken as 1.0 nsec, a temporal 
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coincidence resolution of 4 nsec is achievable with 1 GHz A/D converters and a 1 GHz 
coincidence clock. 

[0057] Knowing the Compton scattering angle determines a cone of origin of the 
incident photon y 0 , as shown in Figure 4A. Once a collection of coincident events has 
been accumulated, a deconvolution is applied to derive the probable source position on 
each cone of origin for each incident photon y 0 , as shown in Figure 4B. Figure 4B 
illustrates multiple coincident events arising from the same point source of incident 
photons, traveling along incident paths 110, 410 and 420, whose cones of origin 433, 434 
and 435 can be deconvolved to the original point source. The paths 114, 414 and 424 of 
the scattered photons are also shown. 

[0058] The distribution of Compton scattered photons from a particular point in space 
can be described as the sum of a set of annula (rings) centered at the intersection of the 
second detection plane with the vector originating at the scattering point and 
continuing along the direction of the incoming photon, as shown in Figures 4C-4E. 
Figure 4C illustrates the probability distribution of absorption events in the second 
detection layer resulting from incident photons traveling along incident path 432 from a 
point source 430 and undergoing Compton scattering at one location 412 in the first 
detection layer. Shown in Figure 4C are the central axis 438 of the forward-projected 
scattering cone from source 430 through location 412, the path 434 of the Compton- 
scattered photon, and the forward-projected Compton scattering pattern 436 from 
source 430 through location 412. 
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[0059] Figure 4D illustrates the combined probability distribution of absorption 
events in the second detection layer resulting from incident photons from two point 
sources 430 and 440 undergoing Compton scattering at the same location 412 in the first 
detection layer. As shown, for the second point source 440 there is a second forward- 
projected Compton scattering pattern 446, defined around central axis 448 from point 
source 440 through location 412. 

[0060] Figure 4E illustrates a skewed probability distribution of absorption events in 
the second detection layer for photons incident at a non-perpendicular angle on the first 
detection layer, and the reconstruction of an unskewed distribution 450 in a virtual 
plane perpendicular to the trajectory of the incident photons (i.e., orthogonal to axis 
448). 

[0061] The Compton scattering angle 6 is a well-known function of incoming 
(incident) photon energy versus outgoing (scattered) photon energy (based on 
conservation of momentum and energy). The distribution of scattering angles as a 
function of energy of incident photons y 0 is described by the Klein-Nishina equation, 
which yields the differential scattering cross-section da/dQ , given by equation (4), 
where a = hv/rn 0 c 2 , and r 0 is the classical electron radius. 



da 



1 



l + a(l-cos0) 



1 + cos 2 e l a 2 (1 - cosO) 7 



(l + cos 2 0)[l + a(l-cos0)] 



(4) 
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[0062] The expected distribution of scattered photons with the same origin, direction, 
and incident energy is, therefore, the set of annula with a radial intensity distribution 
following the Klein-Nishina equation. The measured distribution can, therefore, be 
described as the convolution of the incoming intensity vector and a convolution kernel 
(the scattering kernel), as shown in Figures 6A and 6B. Figure 6A illustrates photons 
along path 610 undergoing Compton scattering in the first detection layer, resulting in 
scattered photons 614 within a range of angles about the original path of incidence 610. 
The forward-projected Compton scattering cone 616 is defined about central axis 618 
from the point 612 of Compton scattering in the first detection layer. 

[0063] Figure 6B illustrates an angular probability distribution of detected Compton- 
scattered photons, in which failure to detect unscattered incident photons or very low 
angle Compton events (where energy absorbed is too small) results in a "hole" around 0 
degrees from the path of incidence 628. The "hole" 624 represents an exclusion zone for 
incident gamma ray photons which did not interact in the first detection layer. Curve 
622 represents the probability of Compton-scattered gamma ray photons absorbed in 
the second detection layer. Confidence limits 626 are optional. 

[0064] Deconvolution is a well-known mathematical method applicable to linear 
distributions. Consider a one-dimensional example of deconvolution in image 
processing, in which an input representing a single point along the axes is "scattered" 
according to a Gaussian distribution. To determine the location of the original point, 
one would replace each value on the line by the sum of each point along the 
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measurement line with the "inverse" of the Gaussian distribution. The inverse can be 
computed or approximated by taking the Fourier transform of the original "scattering" 
distribution, performing multiplicative inversion ( x' 1 ), and then taking the inverse 
Fourier transform, as shown in Figures 6C through 6E. Figures 6C, 6D and 6E illustrate, 
respectively, the inverse deconvolution kernel, the deconvolution kernel, and the 
resulting delta function defining the probable angular origin of incident photons, 

[0065] Another example of a one-dimensional deconvolution is the first step in the 
image reconstruction process employed in transmission computed tomography (CT) 
and emission computed tomography (ECT) imaging. Each one-dimensional projection 
is deconvolved (although it is sometimes referred to as "convolution") prior to being 
back-projected. The back-projection process itself creates the equivalent of a 
convolution of the original projections. Pre-emptive deconvolution allows back- 
projection to result in a correct image. 

[0066] In this particular case, the distribution function is the Compton scattering 
kernel. The Compton scattering kernel is a set of values (for a given input spectrum) 
that has traditionally been computed using Monte Carlo methods. This kernel has been 
used in the forward direction for computing expected radiation dose distribution given 
a specified radiation field, as described by T.R. Mackie et al., "Tomotherapy: a new 
concept for the delivery of dynamic conformal radiotherapy", Med. Phys. 20(6) pp. 1709- 
19 (1993). The inverse of the Compton scattering kernel has been used for 
deconvolution for particular methods of optimizing the specification of radiation fields 
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given a desired radiation dose distribution, as described by T. Holmes et al., "A unified 
approach to the optimization of brachytherapy and external beam dosimetry". Int. J. 
Radiation Oncology Biol Phys. 20(4) pp. 859-873 (1991). The use of a scattering kernel 
which incorporates spectral information has also been described in the radiation 
therapy planning literature, and such a kernel is defined as a polyenergetic kernel by N. 
Papanikolaou et al. "Investigation of the convolution method for polyenergetic 
spectra", Med. Phys. 20(5) pp. 1327-36 (1993). Applying a polyenergetic kernel for 
deconvolution takes into account the probable energy distribution of detected events, 
and allows the CDC technique to require less energy specificity than other Compton 
imaging methods. 

[0067] The Compton deconvolution camera can be treated as a two-step aperture. 
The general formula relating an object and image through an aperture is given by 
equation (5). 

h(r") =f(r)**g(r') (5) 
[0068] In equation (5), f(r) is the source of the emitted photons, g(r') is the aperture 
point spread function, h(r") is the detected image, and each r is the position of an 
incident or scattered photon relative to the line from a source to a given pixel in the first 
detector plane. In the Compton deconvolution camera, the first detector plane acts as 
two apertures. 
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[0069] First aperture: By utilizing coincidence (timing) electronics, the Compton 
deconvolution camera acts like a set of independent pinhole collimators in that the 
second detector plane can be treated as a set of virtual image detection planes, one per 
pixel of the first detector plane. A pinhole collimator has an idealized point spread 
function g(r') = 8(r-r 0 ) . 

[0070] Second aperture: Once the first pixel has detected an incident photon, the 
medium of the first detector plane in that pixel acts as an aperture onto the virtual 
second detector plane with a point spread function defined by the Klein-Nishina 
equation and the likelihood of a priori detection in the first pixel In practical terms, the 
exact form of the distribution in the second detector plane will depend strongly on the 
energy of the incident photons, the construction of the first detector plane (photon 
stopping power determined by atomic composition and thickness), the separation 
between the two detector planes, and the sensitivity of the second detector plane. 

[0071] In addition to the point spread function due to Compton scattering, there is 
also a spatially variant component due to the geometry of the second detector being 
planar instead of a spherical surface centered at the first pixel in the first detector layer. 
This distortion is purely geometric and can be undone (compensated to recover the 
original geometry of the point spread function in a plane orthogonal to the axis of 
Compton scattering) by applying an inverse distortion on the data. The inversion is a 
reverse look-up of D(0,y/), as set forth in equation (6). 
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[0072] In equation (6), w i is the scattering kernel from the Klein-Nishina equation, i 
represents each measured pixel (position) in the second detector plane, y/ is the angle of 
incidence relative to the perpendicular to the two detector planes, and 0 i defines the 
angle between the axis joining the coincidence pixels in the first and second detector 
layers and the axis of incidence. The geometric distortion function can be visually 
identified to be slowly varying monotonic up to the limits of the detection acceptance 
angle (< 80°, but would preferably be further limited to reduce background scattering) 
for a feasible detection geometry. 

[0073] The inverse distortion operation must be applied before deconvolution of the 
image (recovery of the original point spread function) on the second detector plane. 
Recovery from Compton scattering (for straight forward orthogonal incidence to an 
aperture) requires deconvolution in finite terms. A distribution function f(r) may be 
blurred by an aperture function g(r/), resulting in a detected distribution ft(r/') as set 
forth in equation (7). 

Kr") = Yf(n)g(r' l -r l ) (7) 



[0074] Recovery of the original signal f(r) from the detected signal h(r") is 
performed according to equation (8), where G 1 is the deconvolution kernel. 
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i 

[0075] For a linear, shift-invariant aperture g{r'), G' 1 is given by equation (9), where 
the Fourier transform of g(r') is given by equation (10). 

G' 1 = 1 / G (9) 

G = -l=[g{r)e' kr dr (10) 

[0076] Then the inverse Fourier transform of the restoration operation G" ! is given by 
equation (11). 

[0077] The coincidence Compton scattering kernel can be very roughly approximated 
for demonstration purposes to be g(r') as set forth in equation (12), where a is the 
variance of Compton scattering angles. 

g( /) = \/\ e -^ 1 (12) 

[0078] This is a Gaussian distribution with a hole in the center (Figure 6D) due to 
failure to register a coincidence event when the energy deposited by the incident 
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photon in the first detector layer is below the threshold of detection (includes the case of 
no Compton scattering at all). 

[0079] Given C, the contraction operation, g(r'), the scattering kernel, and /z(r".), the 
detected signal (Figure 6C), for pixel i in the first detector plane, the recovered 

sinogram (projection) image fir) (Figure 6E) on the virtual detector plane is as given in 
equation (13). 



[0080] This sinogram image is then used to determine what values to accumulate in 
the single-photon emission computed tomography (SPECT) sinogram space (a set of 
planar sinograms). 

[0081] The detection geometry defines the relationship. For a detector at some 
nominal angle with respect to the vertical (e.g., the gantry angle), the incoming angle of 
the detected photon streams is calculated from the relative positions of the first pixel in 
the first detector plane for each recovered sinogram (projection) pixel in the second 
detector plane. The angular offset iff is the equivalent rotation of the detector heads for 
sinogram construction. From this, the offset x from the center of an equivalent parallel 
collimator camera can be computed. 

[0082] Figure 6F illustrates the first and second detector layers (662 and 664) rotatable 
about an isocenter 660. A point radiation source 670 emits a y 0 photon incident along 




(13) 
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path 672. This photon is detected in the first detector layer at a displacement y ' relative 
to the detector radius vector 674. The axis of incidence intersects the second detector 
layer at a displacement y ". The first and second detector layers (666 and 668) are 
shown in a second position with radius vector 676, chosen for rebinning purposes to be 
perpendicular to this incident path (for the purpose of equivalent parallel coUimation). 
In this second detector position, the same incident y 0 photon would be detected in both 
detector planes at a new displacement x. 

[0083] The angle ^and offset x values for these parallel collimator equivalents in two 
dimensions are set forth in equations (14) and (15), respectively. 



[0084] In the above equations, d n is the separation between the two detector layers; R Q 
is the radius from the isocenter to the first detector layer; y/ is the angle of incidence 
relative to the detector radius vector; <j> is the angular displacement of the pixel in the 




(14) 



first detector layer relative to the detector radius vector when it is in the first position; 



and equation (16) defines the relationship between y', R 0 and 0. 



tan" 1 — 
Ro 



(16) 



27 



[0085] The sets of projections can be separated (or associated) by angle of incidence of 
the incoming photons (rebinned). The incoming photon angles will already have been 
estimated and tabulated as a result of the deconvolution of the measured scattered 
radiation field by the Compton kernel. While a (physically) collimated SPECT camera 
is restricted to measuring a single set of projections for a given position of the 
detector (s), the CDC simultaneously measures multiple projection sets. 
[0086] Having rebinned all of the data into a set of sinograms (projections) 
corresponding to a parallel collimated SPECT scanner, the image is reconstructed using 
a standard SPECT image reconstruction algorithm. Judicious choice of binning 
geometry also allows for use of Fourier reconstruction techniques. 

[0087] The volumetric placement of counts provides a three-dimensional 
reconstruction that can be used for planar projection displays, or may be used in 
combination with rotational SPECT to provide a more accurate three-dimensional 
dataset. Because the incident photons y 0 are arriving from different directions, even 
planar image reconstruction requires estimation of depth of origin. This is also 
accomplished by three-dimensional reconstruction based on all vectors within the field 
of-view defined by the acceptance angle. 

[0088] Because the two-dimensional angle of origin of each the incident photon y 0 is 
known, coincident counts can be selectively accepted according to angle of origin to 
generate collimated, steered, focused, and sector images. This provides substantial 
flexibility in image reconstruction and reprocessing in software post-processing. 
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[0089] Figure 8 illustrates the overall process of performing event detection and 
image reconstruction using the CDC, according to one embodiment. Ionization (or 
scintillation) events are detected in the first and second detection layers at blocks 810 
and 820, respectively. The detected events are pre-amplified and A/D converted for the 
first and second detection layers at blocks 812 and 822, respectively. At blocks 814 and 
824, the two-dimensional positions (in the plane of the detection layer) of the events are 
determined in the first and second layers, respectively, using two-dimensional position 
interpolation. 

[0090] Next, at block 830 a determination is made of whether there are two or more 
events in coincidence. If not, the process loops back to blocks 810 and 820. If there are 
exactly two events in coincidence, the process continues from blocks 816 and 826, in 
which the energy e re of the recoil electron for the event in the first and second layer, 
respectively, is estimated, using spatial integration of pulse height. If there were more 
than two events found to be in coincidence at block 830, then the process continues 
from blocks 816 and 826 only if there were exactly two events found in coincidence in a 
single zone. Otherwise, the process loops back to blocks 810 and 820. 

[0091] After estimation of the recoil electron energies e re , the estimated energy values 
£ re are applied to a gating function at block 834 to validate that the two coincident events 
are in fact related to the same incident gamma photon. If the validation fails, the 
process loops back to blocks 810 and 820. If the energy values are validated, then at 
block 836 a determination is made of whether the energy e o of the incident photon 
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equals the sum of the energy e re of the recoil electron and the energy of the scattered 
photon. If such equality is found, the process continues from block 838. Otherwise, the 
process loops back to blocks 810 and 820. 

[0092] At block 838, an appropriate number of counts is stored for the detected 
coincident event in each detection layer in the proper location (pixel) of memory for 
each layer. Next, at block 840 the Compton scattering angle is computed from the 
energy e x of the scattered photon in the second detection layer. If sufficient counts 
accumulated according to a predefined count threshold or acquisition time limit at 
block 842, then the process continues from block 844. Otherwise, the process loops back 
to blocks 810 and 820. At block 844, the cones of origin are deconvolved into probable 
source locations for the incident photons. If collimation is selected at block 846, then 
any skewed vectors are eliminated at block 850, and a two-dimensional planar image is 
displayed at 852. If collimation is not selected, then a tomographic three-dimensional 
image is reconstructed at 854, and two-dimensional "slice" images or volume 
projections are displayed at 856. 

[0093] Thus, a method and apparatus for generating images in a medical imaging 
system based on detected radiation by using Compton scattering principles have been 
described. Although the present invention has been described with reference to specific 
exemplary embodiments, it will be evident that various modifications and changes may 
be made to these embodiments without departing from the broader spirit and scope of 



30 



the invention as set forth in the claims. Accordingly, the specification and drawings are 
to be regarded in an illustrative sense rather than a restrictive sense. 
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